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Abstract. Inwardly curved polymer brushes are present in cylindrical and spherical micelles or in mem- 
branes tubes and vesicles decorated with anchored polymers, and influence their stability. We consider 
such polymer brushes in good solvent and show that previous works, based on a self-similar concentric 
structure of the brush, do not describe the most stable structure. We use scaling laws to derive very simply 
the leading term of the free energy in the high curvature limit, where the osmotic pressure is the relevant 
physical ingredient. We also derive the complete conformation at all curvatures using a self-consistent field 
approach. The free energy is computed therefrom using a local scaling description 

PACS. 3 6.20.Ey Conformations (statistics and dynamics) - 82.35.Gh Polymers on surfaces, adhesion - 
82.70.-y Disperse systems; complex fluids 



1 Introduction. 

Various types of polymer layers have been studied in re- 
cent years. The most common situations include polymer 
brushes (where polymers molecules are end-tethered to a 
repulsive surface or interface), reversible adsorption from 
a dilute polymer solution, irreversible adsorption from a 
melt. These are important in the field of adhesion, wet- 
ting, lubrication, colloid stabilization Different the- 
ories have been proposed to describe the structure and 
' the properties of the polymer coatings as a function of the 
, mode of attachment, the quality of the solvent, etc. How- 
ever, concerning the geometry of the surface which sup- 
ports the layer, most of the work is concerned with flat or 
convex interfaces This is in contrast with a number 
of situations where polymers chains are attached to con- 
cave surfaces. For example, diblock copolymers can form 
spherical micelles in a selective solvent. The core of these 
is made of the segregated moeties and can be pictured as 
a dense layer of chains "grafted" to the inner surface of a 
sphere. Recently, nine different phases have been found in 
a ternary system of poly (ethylene oxide)/poly(propylene 
oxide) amphiphilic block copolymer, water and oil (selec- 
tive solvents). Six of these phases are made of cylindrical 
and spherical normal micelles (water-in-oil) or reverse mi- 
celles (oil-in- water) ||. Lattice Monte Carlo simulations 
were performed to study chain conformations in spher- 
ical cavities 0. Polymers attached to fluid membranes 
provide another interesting situation. Indeed, in biological 
systems, lipid bilayers are often decorated by a large num- 
ber of macromolecules. Inspired by biological structures, 



simplified model systems of polymers attached to lipid bi- 
layers have been studied experimentally pp| JlO| , pd[ | . It was 
shown that anchoring polymers induces dramatic changes 
in the vesicles || or multilamellar membrane tubes [ pd] 
shapes. In general, the polymers will be attached to both 
sides, i.e. to both monolayers of the lipid or surfactant 
bilayer. Anchored polymers also stabilize the equilibrium 
formation of small unilamellar surfactant vesicles |Q over 
a wide range of composition and temperature range: the 
vesicles are spontaneously spherical and monodisperse and 
their radius is small (comparable to the natural size of the 
anchored polymers). This spontaneous tendency for bilay- 
ers to bend so to form vesicles has been explained by a 
coupling JT^ I between the curvature and the asymmetry 
of grafting polymers on each monolayer fl4|| . 

Polymer adsorption finds extensive applications in mod- 
ifying physical and chemical properties of surfaces. In this 
context, experiments have highlighted the effect of surface 
roughness at the size scale of the grafting polymers on the 
total amount of polymer adsorbed [ fl5| . A significant en- 
hancement of adsorption was observed when the bead to 
polymer radius is of order of unity, indicating that poly- 
mer brushes will prefer to adsorb on spherical or cylindri- 
cal bumps (i.e. defects of large positive curvature), since 
the deformation energy of the brush will be lowered in 
these geometries with respect to the flat case @. On the 
contrary, valleys, (i.e. defects with large negative curva- 
ture) substantially reduce the grafting. Such a reduction 
was understood in terms of topologically reduced steric 
hindrance to adsorption. Monte-Carlo computer simula- 
tions were also performed |l7| to study the adsorption of 
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triblock copolymers on rough surfaces and confirm that 
grafting on defects of negative curvature (concave sur- 
faces) is not favoured because of conformational energy 
penalty of the anchored polymers. 

Surprisingly, even if it was anticipated that in the high 
curvature limit, concave and convex brushes should not be 
described in a similar way jl6| , very few theoretical works 
have been devoted to the specific case of inwardly curved 
brushes up to date. To our knowledge, only a few theoreti- 
cal papers have addressed this problem, two of them in the 
context of the study of the elastic properties of polymer- 
decorated membranes 0,^9) , one in the context of shear 
swelling of polymer brushes grafted onto convex and con- 
cave surfaces pCfl , two in the context of microemulsion 
stabilization by diblock copolymers |2l| , p2| . E. Zhulina p3| 
also pointed out references P4p3] . However, most of these 
authors described the chain conformation and the free en- 
ergy of the brush on the concave side of the interface by 
using an inverted version of the classical Daoud-Cotton 
scaling approach pf|. In this approach, the structure of 
the brush is self-similar and scales linearly with the dis- 
tance from the center of curvature (Figure 1). It correctly 
describes brushes anchored on the convex side of interfaces 
but not on concave side, because, as we will show below, 
this structure is not the most stable one. 

The main concern of the present work is to provide a 
correct description of the concave brush structure. In or- 
der to clarify the situation of interest, let us briefly review 
the different conformation regimes (Figure 2), depending 
on the number / of polymer chains (assumed to be linear, 
flexible, made of identical, neutral monomers of size a), on 
the number N of monomers per chain and on the radius 
R of the inner surface of the sphere. The number of chains 
per unit surface area of chains (grafting density) can be 
written adimensionally as a = a 2 f /(AttR 2 ). For simplicity 
reasons, we only treat the case of an athermal solvent : the 
excluded volume parameter is v = 1 (in units of a 3 ). We 
assume that the chains are repelled by the interface. An 
important practical issue requiring attention is the case 
of flexible membranes with grafted polymers. From this 
point of view, we will now express all quantities in terms 
of the grafting density, a, rather than in terms of the num- 
ber of chains, /. It allows also an easier comparison with 
the limit case of a flat brush. 

At large sphere radius R, if chains are grafted too far 
apart to interact, they form separate mushrooms (regime 1 
of Figure 2) whose Flory radius of gyration is Rf = aN 3 ^ 5 . 
They start to interact when a = (a/Rp) 2 , i.e., when : 

a ~ N~ 6/5 (1) 

For larger values of the grafting density a, chains are 
stretched in good solvent and form a brush whose height 
is essentially the same as in the planar geometry : 

Lfiat^aNa 1 ' 3 (2) 

This is regime 2 of Figure 2. It is limited by two phenom- 
ena. First, the grafting density cannot exceed one chain 



per squared monomer size, for obvious steric reasons. This 
condition reads : 

o- < 1 (3) 

Secondly, the sphere curvature becomes important when 
R ~ Lfi a t, i.e., when : 

R ~ aNa 1/3 (4) 

At smaller radii, the description of the grafted chain con- 
formation is not straightforward since both the curvature 
of the sphere and the interaction between neighbouring 
chains are important. These are regime 3 (where chains 
are still stretched) and regime 4 of Figure 2 (where chains 
are confined). The separation of these two regimes occurs 
when the radius is equal to the natural size of a chain in 
a semi-dilute solution (see Section 3.2), i.e. when 

R ~ aa-^N 3 ' 7 (5) 

They are limited by two conditions: one occurs when no 
more solvent is present : 

R > aNa (6) 

and the second one when there is only one chain in the 
sphere : 

R > aa- 1 ' 2 (7) 

In regimes 2 and 3, chains are stretched and their con- 
formations should play an important role in the struc- 
ture and the energy of brushes in these regimes. We thus 
briefly review two usual types of descriptions. In scaling 
descriptions, all chain conformations are identical and in 
particular, chain ends are all located at the same distance 
from the grafting interface. Conversely, self-consistent field 
(SCF) descriptions allow for a spatial distribution of chain 
ends. Correlatively, all chain conformations are not iden- 
tical : at any given distance from the wall, chains are 
stretched differently according to how far from the graft- 
ing surface their free ends are located. 

The common feature of both approaches is that the 
brush conformation results from the interplay between 
the two-body interactions between monomers pertaining 
to the brush chains (which tend to stretch the brush fur- 
ther out) and the elastic stretching of the brush chains 
(which limits the brush extension). As long as these two 
components are well taken into account, both types of ap- 
proaches yield essentially similar conclusions. 
Thus, the Daoud-Cotton description of convex brushes, al- 
though based on scaling arguments, yields correct results. 
Let us emphasize a few features of such a system. In real- 
ity, free chain ends can sustain no stretching, whereas the 
grafted ends are under tension. In the model by Daoud 
and Cotton, the typical blob size (see Figure la) is larger 
at greater distances from the wall (it is proportional to the 
distance from the center of curvature). Thus local stretch- 
ing is weaker far from the wall. This feature ensures local 
equilibrium of the chain strands: the gradient in concen- 
tration and thus in two-body interactions tends to push 
them away from the wall, whereas the gradient in degree 
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of stretching tends to pull them back towards the wall. If 
we focus on the last blob, the elastic force, given by the 
Pincus law ]27[| , is ksT/^(L), where £(£) is the size of the 
last blob. Since it is large, the elastic force is small and is 
balanced by the small concentration gradient at the edge 
of the brush. The success of this scaling vision is thus due 
to the fact that chain-ends are segregated in the last blob 
and they sustain almost no stretching. Hence, the basic 
ingredients are thus explicitely present in this scaling de- 
scription. 

Convexity of the interface is a key feature in this success 
of the Daoud-Cotton model : the planar geometry is a 
marginal case for the Alexander-de Gennes scaling model. 
It yields essentially the same results (brush extension and 
free en ergy ) as the pionneering SCF description by Se- 
menov p8j which was further developed independently by 
two groups ||j||3(J. However, the monomer concentration 
profile found in the SCF theory is parabolic and vanishes 
at the end of the brush. This is different from the step 
profile of the scaling description which imposes a sharp 
density gradient at the edge. We might argue that this 
gradient balances the elastic restoring force exerted on 
the stretched chain and that it ensures the global equilib- 
rium of chains. But it has been shown that this is not the 
most stable structure since the SCF theory yields a lower 
minimum free energy. 

Things get worse in the case of a concave geometry. 
Two scaling descriptions can be found in literature : one, 
adapted from the model by Daoud and Cotton |l8[pj| , pi| , p4[ , 
the other, given by Sevick [|o|, is based on an Alexander- 
de Gennes brush treated locally to account for a cylindri- 
cal geometry. The essential feature of the "inverse" Daoud- 
Cotton approach is that the blob size is still taken as pro- 
portional to the distance from the center of curvature : 
because the interface is concave, the blob size decreases 
as we go further away from the wall (Figure lb). The es- 
sential results of this model are as follows (Throughout 
this article, free energies F are in units of thermal energy 
k B T) : 

— the thickness of the brush is L DC ^ R[l-(l- aN ° 1/a ) 3/5 } 

— the free energy per chain is Fbc — —Act 5 / 6 ^^j^r/STR 

2/3 

— the volume fraction profile is cf>(r) ~ nrzTTRWa wnere 
r is the radial distance to the spherical surface (see 
Figure Ibis) 

Among these predictions, the last one is certainly the 
most surprising. It predicts that the volume fraction in- 
creases when we move off from the surface, and notably, 

2/3 

at the edge of the brush, we find 4>{L) ~ w^^xjwjmiji ■ 
Thus forces exerted on the last blob (osmotic gradient and 
elastic force) are very important. Eventually, the volume 
fraction diverges at the center for R ~ a Act 1 / 3 . Hence, 
even if the local stability inside the brush is fulfilled, as 
in the convex case, we might expect that this profile does 
not give the best answer. 

Among the results yielded by the Sevick model, one 
is particularly surprising: at constant grafting density, the 



brush height is predicted to decrease if the interface is 
curved more strongly. Yet simple arguments indicate that 
on the contrary, the brush height should increase. Suppose 
for a moment that the height does not change when the 
curvature is changed. Curving the interface more strongly 
then reduces the overall volume of the brush. Two-body 
interactions are then more important, which tends to swell 
the brush. The height will therefore stabilize at a some- 
what greater value. 

A third approach considers the droplet curvature effect 
on the emulsion stabilization by diblock copolymers p2[ . 
They study the small corrections due to curvature to the 
flat case in the framework of the Alexander-de Gennes 
model, which is valid for very weak curvatures. 

To describe such concave brushes under large curva- 
tures, one must therefore keep in mind that the curved 
brush is some kind of a compressed brush. Curving the 
interface indeed amounts to increasing the brush concen- 
tration : the corresponding compression does not originate 
in an interaction at the brush edge, however, but rather 
in some kind of lateral interactions (and the brush profile 
is very different from that of a compressed planar brush, 
as we indicate in Appendix A). In this context, one result 
obtained by Rabin and Alexander on compressed polymer 
brushes is most important (3^]: they showed that, when 
the layer thickness is forced to be smaller than that of a 
free brush (which here corresponds to regimes 3 and 4), 
the main contribution to the free energy of a compressed 
brush in good solvent is the osmotic contribution. The 
elastic stretching energy of the chains is comparatively 
smaller, even though they are still strongly stretched (at 
least in regime 3). As a result, in the limit of strong com- 
pressions (regimes 3 and 4, R <C ciAct 1 / 3 ), the osmotic 
pressure is the mostly contributing term and is therefore 
uniform over the whole brush volume. In other words, the 
monomer concentration is uniform over the entire brush 
and the brush structure is locally that of a semi-dilute 
solution of equal overall concentration. 

The size of the semi-dilute blobs is £ ~ a</> -3 / 4 (where 
4> — 3aNa/R is the volume fraction) and the free energy 
of the whole brush is therefore : 

F 3 s P h 4 ~|J~(ACT) 9 / 4 i? 3 / 4 (8) 

This expression for the energy is valid well inside regime 3 
or in regime 4, i.e., when the brush is strongly compressed 
(R^L f i a t). 

Our aim is to give a correct description of the chain 
structure and free energy in regimes 2 and 3. The present 
work is organized as follows. 

We present in Section 2 a SCF approach for such con- 
cave brushes valid for a near-theta solvent and an alterna- 
tive convolution method is outlined in Appendix A. The 
purpose of Section 3 is to give a scaling picture of concave 
brushes by refining SCF results taking into account cor- 
relations between monomers. A discussion of lateral fluc- 
tuations of grafted chains is given in Appendix B, where 
it is shown that they do not alter significantly the chain 
free energy and thus do not endanger the SCF approach. 
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In Appendix C, we give the results obtained in the cylin- 
drical geometry through a scaling approach. 



profile is parabolic and we know the form of the local 
extension of the chain : 



2 Self-consistent mean field description. 

In regimes 2 and 3, the polymer chains are strongly stretched 
and it is thus possible to use the self-consistent field method 
p8| , p9|j3C| ] . We note r the radial distance to the spherical 
surface[defined by r = 0). 

The aim of this section is to determine the free en- 
ergy F by solving for the monomer volume fraction <p(r) 
and for the distribution g(r) of chain-end monomers self- 
consistently. In the classical limit of strong stretching 
the free energy of the system is 



F = 



R pro 

9(ro) / 
o Jo 

1 r R 4w(R 



o 



a- 



3 
2a 



") 



e(r, ro) 



drdrn 



(9) 



where g(ro)dro is the number of chain ends in the spherical 
shell of radius R — tq and width dro and e(r, ro) = gjj 
is the local extension of a chain whose free end is at ro 
(the variable n denotes the current monomer). The first 
term corresponds to a sum of single-chain free energies in 
the potential fcsTw0(r), the second term corrects the fact 
that monomer-monomer interactions have been double- 
counted. The normalization conditions 



dr 



ro 

(i e(r,r ) 



= N 



(10) 



l (r)=<P 



G(r- xR) 



TT 2 R 3 



where 



B 



R* 



10 / R 



8aN 3 a 3 v 9 \ R* 

1/3 



9tt ; 



a(va) 1/3 N 



(13) 
(14) 

(15) 
(16) 



where is a step function, $ = 3aNa/R, B is obtained 
from the equal-time argument (Eq. (|To|) ) and A is related 
to x by (Eq. ®) 

(3 - 3* + x 2 )xA = 1 + y (§) 3 * 3 (1 ~ \x + \x 2 ) (17) 



The minimization of the free energy per chain F(x) = 
F[(f> m i n ,x]/ f with respect to changes in x under the three 
conditions a) < x < 1, b) (j>(x~) > and c) 0(0) = 
(PA < 1 is illustrated in Figure 3. This gives rise to a 
distinction between two cases depending on whether R is 
smaller or larger than R* , where R* (see Eq. (|l4|)) is the 
critical radius for which the sphere is entirely filled up by 
the brush: 

— for R > R* (regime 2), the minimum free energy F 
occurs for x* defined has the unique solution of the 
equation 



4ir(R - r) 2 (j)(r)dr = 4irR 2 aNa 



(11) 



fix the number of monomers in one chain and in the whole 
system. The mean field concentration is determined by 
5F/Scj) = 0, which yields 



4n{R 



g( r o) 

e(r,r ) 



dr 



(12) 



We perform this calculation by two steps [£9|: first 
we suppose that the brush height, xR, is fixed, and we 
minimize the free energy given by Eq. (|TJ) with respect 
to changes in the field </>. This yields a family of profiles 
<j>min- Then we determine this unknown parameter x (with 
< x < 1) by minimization of F[(p min , x] which yields the 
brush height, x* , at equilibrium. 

The "equal-time" property can be applied for the self- 
consistent field as was done previously in the case of a 
melt |28| if the following two assumptions hold : a) there 
is no dead zone containing no chain ends fl32"| , which can 
be written as g(r ) ^ for < Tq/R < x; and b) the chain 
ends are not under tension : e(ro, ro) = 0. From the equal- 
time argument, we know that the monomer concentration 



3 I 4 2 5 
x 3 [ -x x 



9 



20 



(18) 



such as condition a) is fulfilled (Note that this condi- 
tion is equivalent to the smooth vanishing of the poly- 
mer concentration at the edge of the brush [Q). It 
remains a central region of radius R(l — x*) in the 
center of the sphere which contains pure solvent ; 
— for R < R* (regime 3) the minimum of F occurs for 
x = 1, there is no empty region and the concentration 
at the center of the sphere is finite. 

By varying R from oo to aN 1 / 2 , the limit of validity of 
the SCF approach, the following scenario emerges. 



2.1 Regime 2 (R > R*) 

In regime 2, the brush height x* is given by Eq. ( |l8| ) and 
the monomer volume fraction (Eq. (|l3|)) simplifies into 
H (Figure 4) 



(r) = <PB 



(19) 
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Equations (112) and (fL9|) yield the chain end distribution 

— — ( i 9 

written in units of chains per unit length g(t) — SLM 



AirRa 



3/2 



g(t) = Bx*H [3(2 - 3x* + 2x* 2 )VT~^t 2 - 8x* 2 {l - t 2 ) 
+6x* arg tanh a/1 - t 2 - 9x*t 2 arg tanh yjl-t 2 \ 

(20) 

where t = ro/(x*R). The distribution of chain ends can 
also be expressed in terms of a number of chain ends per 
unit volume : p{t) = jsnz^rp (in units of chains per unit 
volume, <P/N) for < t < 1. A typical probability density 
for free chain ends is shown in Figure 5. 

^From these, we readily derive the free energy per 
chain F 2 = F(x*) given by Eq. (Jsl) : 



which is also the classical result. 

Going further, we can expand this result for very small 
curvatures. Equation ( |l8| ) yields the asymptotic value of 
the brush height 



J flat 



1 + T 



1 Lfi a t 

4 R 



29 
240 



-'flat 

~R~ 



(25) 



Note that a simple scaling method does not yield the 
correct value for this asymptotic increase in the brush 
height |3?J . We can also expand the free energy for small 
curvatures 



F 2 _ p. 



F 2 = 



_9_ 

10 



1/3 



2flat 



N{vcr) 2/3 lP 



R* 



5 L 



flat 



12 R 



73 
336 



-'flat 

~R~ 



(26) 



(21) 



where 



10. 4 1 / 3 



3(3 - 3x* 4 
(25/6) 1 / 3 x 
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1/3 



(25/6) 
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201 
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We observe that this is consistent with the results of Mil- 
ner and Witten |38| who found for the convex geometry 
the same development by replacing R by — R. This is be- 
cause, as explained in |38|] , in the convex case, the dead 
(2§tpe yields negligible terms for small curvatures p9| |. 

In conclusion, for this region 2, the SCF calculations 
show a weak increase of the brush height when the inter- 
face is curved more and more, which involves an increase 
of both elastic and osmotic energy. 

*4\ 



In the entire regime 2, equation ([18]) shows that the nor- 
malized brush height x* increases when the curvature R* / R 
increases, until it reaches x* = 1 for R — R* (Figure 6). In 
other words, the brush height increases when the interface 
is curved more strongly, as long as there is some free sol- 
vent in the center. In the same way, it can be shown that 
the normalized free energy tp (Eq. (p3|)) increases mono- 
tonically at constant grafting density when the curvature 
increases (Figure 7). 

The osmotic pressure balances the elastic restoring force 
since the polymer brush is at equilibrium under no exter- 
nal compression. This readily explains HTJ] the fact that 
the monomer volume fraction vanishes smoothly at the 
outer edge of the brush (equation ( p"8| ) is equivalent to 
c/>(x*)=0). 

Izzo and Marques Q seem to have also used an SCF 
approach for swollen concave brushes in this regime. They 
find the same equation as Eq.(|l8|) for the brush height, but 
their free energy is different from the one given here. We 
do not understand why at this stage. 

In the limit of zero curvature, i.e. R — ► oo, the brush 
height L 



2.2 Regime 3 (aN 1 / 2 < R < R*) 

In regime 3, the brush height L is fixed at R (x = 1). The 



monomer volume fraction becomes |35j (Figure 4) 



(r) = cj){R) + <PB 



1- 



,r 
\R 



(27) 



where 'P(R) is the volume fraction 
sphere given by : 



4>{R) = $ 



1-1^ 

R* 



at the center of the 



(28) 



which increases when R decreases. Thus, the monomer vol- 
ume fraction tends to become more uniform as the radius 
R is decreased. 

The chain end distribution is (t = r/R) 



g(t) = t [-6 ( 
+6(1 



1 



^/l^t 2 -8B(l-t 



2)3/2 



x*R given by Eq 

1/3 



Lflat = 



R* = 



(18) can be expressed as : 

1/3 



7T Z 



a(va) 1/3 N (23) 



jqB) arg tanh Vl - t 2 - 9Bt 2 arg tanh Vl - t 2 \ 

(29) 

As is shown in Figure 5, the major number of chain ends 
per unit surface are located at the edge of the brush. The 
free energy per chain is given by : 



which is the result of Milner, Witten and Cates 
over, we have ?p(0) = 1 and the energy Eq. (|2l| 



F- 



2flat 



1/3 
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R_ 

R* 
(30) 

Mathematically, this free energy F3 is frozen at F(x = 1) 



because the equilibrium state x 



> 1 is out of reach, 
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x = 1 being a physical bound of the problem. This situa- 
tion is similar to the case of flat brush under compression 
(a pressure is exerted at its edge) : the concentration pro- 
file has a truncated parabolic shape and its value at the 
edge of the brush, Eq. (p8|), increases when the compres- 
sion is stronger [^0). As it has been shown by Milner et 
al., this means that there can be no "interdigitation" of 
chains at r = R |34|j . 

At constant grafting density, F 3 (Eq. (|(j|)) is a strongly 
increasing function of R* / R (Figure 7). A careful analy- 
sis of this energy shows that for R <C R* , the last term 
is negligible, it comes from the parabolic nature of the 
monomer concentration. In the same way, from Eq. ( p8| ) 
we can reasonably consider that the monomer concentra- 
tion is constant and equal to $ = 3Naa/R. Therefore, 
the free energy can be seen at the level of a simple Flory 
theory 

~ aN 2 a R 2 

F 3 ~ v—pr- + (31) 



R 



Na 2 



where the first term is the osmotic pressure (at a mean- 
field level) and the second one the elasticity of an ideal 
chain. For R = R*, these two energies balance, but when R 
decreases further, the excluded volume interactions strongly 
dominate and the energy is practically equal to the os- 
motic pressure (Figure 7). The point is that these two 
energies do not balance because the brush height is geo- 
metrically constrained. 



3 Scaling refinements 

In the previous section, we used a self-consistent field ap- 
proach to study the conformation of a concave brush. Such 
an approach describes both the average distribution of 
free chain ends and the average position of any monomer, 
conditionally if the position of the end-monomer of the 
chain to which it pertains is known. It is known classi- 
cally in two versions. The mean-field version such as that 
exposed in the previous section provides a somewhat er- 
roneous expression for the brush free energy, which varies 
like 4> 2 . This flaw is intrinsic to the mean- field approxima- 
tion and is present also in the Flory approach. A correct 
expression for the free energy can be obtained using scal- 
ing concepts JIT) and varies like tfi 9 / 4 . It was first used 
by de Gennes in a scaling SCF approach for the problem 
of polymer adsorption J42{ . Milner, Witten and Cates ap- 
plied it in order to work out the conformation of a flat, 
swollen polymer brush in good solvent H]. The result- 
ing brush concentration profile obtained after minimiza- 
tion is slightly different from that obtained through the 
mean-field SCF calculation. This is in contrast with the 
original Alexander-de Gennes model for the brush where 
both the mean-field and the scaling forms of the free en- 
ergy yield the very same brush conformation, whose con- 
centration profile is uniform. In the SCF approach, the 
non-uniformity of the concentration is the origin for the 
slight difference between mean-field and scaling flis)] . In 
this section, we present these scaling tools for the vari- 
ous regimes of the spherical concave brush. For stretched 



chains (regimes 2 and 3), we simply calculate the scaling 
free energy from the brush conformation determined from 
the mean-field SCF approach exposed in the previous sec- 
tion. Although it is not as rigorous as the calculation of 
reference ]34|], it should be a sufficient approximation at 
least in regime 3 where the concentration profile is essen- 
tially uniform, as we just mentioned. In any case, the value 
of the free energy obtained in this way in regimes 2 and 3 
will be improved as compared to the mean-field value ob- 
tained in the previous section. 

We also describe more fully chain fluctuations around 
their average conformation given by the SCF treatment. 
Confined chains (regime 4) cannot be treated by the SCF 
theory : we describe the chain conformations and calculate 
the scaling free energy. 



3.1 Scaling refinements in regimes 2 and 3 

The scaling form of the free energy given in reference [ [34] 
is : 



o pR pr 

+ 2a 2 ~J dr °- 9(ro) / < r ^o)^ 1/4 (r)dr (32) 

where the factor 9 / 4 in the osmotic term comes from the 
ksT per blob Ansatz, while the (ft 1 / 4 factor in the elas- 
tic term originates in the correction to the unperturbed 
conformation of the chains which are swollen below the 
length scale £(r) = a<^~ 3 / 4 . The authors interpreted the 
elastic term as a perturbation of each blob of size £ : it 
corresponds to kT (Z\£/£) 2 per blob, where A£ is the aver- 
age blob elongation. But this term can be interpreted in a 
more intuitive way, using a more precise description of the 
local chain conformation. As mentioned briefly by the au- 
thors, at small length scales, chains are like in a semi-dilute 
solution (swollen below £ and Gaussian above £), while 
they are stretched at large length scales. In other words, 
the root mean square distance between two monomers n 
and n' along the chain (see Figure 8a) is given by : 

y/( r _ r ')2 („_ n ')3/5 (y / ( r _ r /)2 <£ = a 0-3/4) 

oc(n-n') 1/2 (f < ^J^Tf < A) 

(y '{r - r') 2 > A) (33) 



oc n — n 



where the intermediate length scale A, where this transi- 
tion occurs, is called the elastic blob (see references 44 4f|) 



A(r, r ) =i a 2 



dr 
dn 



(.-1/4 



(34) 



At this length scale, the polymer chain is not an ideal 
random walk any more : instead, it has no choice but to 
move forward at every step (exponent 1). Hence, the chain 
is mainly a linear string of these elastic blobs stretched 



M. Manghi et al.: Inwardly Curved Polymer Brushes : Concave is not Like Convex 



7 



along the radial direction; the important point being that 
these blobs arc defined along one direction and that in the 
other two directions, blobs of different chains overlap. 

In the elastic term of equation @ , the degree of stretch- 
ing e(r, r*o) must be replaced by an effective degree of 
stretching e e g(r, ro) = a 2 / A(r,ro), which yields the fac- 
tor cj) 1 / 4 in equation (|||). 

Note that the osmotic blob size £(r) depends only on the 
position r, whereas the size of an elastic blob A(r, ro) 
on a particular chain also depends on the local degree 
of stretching of that chain, i.e., on the position ro of its 
free end. 

In principle, the chain conformations should be ob- 
tained through minimization of equation ( 3^) , as in refer- 
ence |34| . Here, however, we simply use the volume frac- 
tion (f>(r) (equation ( p~s| ) or (p7|)), the free end distribution 
g(t) (equation ( po|) or (29)), and the degree of stretching 
derived through the mean-field method developed in the 
previous section (Eq. |l6| ). 

The free energy of the brush in regimes 2 and 3 can be 
computed numerically from equation ( |32"|) along the lines 
we have indicated. 

Here, for simplicity, we focus on the limit where R <C 
R* (well in regime 3). The free end distribution given by 
equation ( p9[ ) at zero order in R/R* then becomes : 

g(t) ~ 6i jargtanhVl - t 2 - yjl -t 2 } (35) 

In the same way, at this order, cf>(r) is constant where 
4>{r) = <P = ZaNa/R. With these two approximations, 
the total free energy is equal to the first two terms of Eq. 



3.2 Confined chains (regime 4) 

The natural size of an iV-monomer chain in a semi-dilute 
solution of volume fraction <j) is known to be R{4>) — 
aAr i/2^-l/8 = £(N/ g )i/* where £ = a ff 3 / 5 = a^" 3 / 4 is 
the osmotic blob size defined earlier in this section. In 
regimes 2 and 3, when the sphere radius is decreased, the 
chains remain stretched as long as their end-to-end ex- 
tension is larger than the elastic blob size A. This breaks 
down when the radius of the sphere reaches the natural 
size of the chain, i.e., R — R(4>) or : 



R 



-l/7 N 3/7 



(37) 



For smaller sphere radii (which corresponds to regime 4, 
see Figure 3), the grafted chains are confined. More pre- 
cisely (see Figure 8b), the chains are swollen at small 
length scales, Gaussian at length scales between £ and 
R, and confined (zero exponent) for larger numbers of 
monomers |Q. The transition between the Gaussian and 
the confined statistics corresponds to G — § X I 4 R 2 /a 2 mono 
At scales larger than G monomers, because of the spherical 
boundary, the chain has no choice but to stay in the same 
region of space instead of performing a random walk : 
globally, it is a string of N/G overlapping blobs of size 
R. Using the same argument as earlier in this section, we 
conclude that the compression contribution to the overall 
energy can be estimated as : 



4el 



N 
G 



N^ 4 a^ 4 (R/a)- 7 ^ 



(38) 



( ^0| ) multiplied by 



F ~ - 



3 aN 2 va 
2 



,1/4 



.1/4 



3tt 2 R 2 



This, however, is only a correction to the free energy, 
whose main contribution is osmotic and is given by equa- 
tion (§). 



a/4 



R 

1/4 



9tt 2 
10 



80 

5/12 



No 2 



(36) 



— a 5 / 6 iV 



5/4 



R 

R* 



7/4' 



where the dominant contribution is osmotic (first term). 
The second term is the elastic term (computed above for 
R/R* — > 0). The next nonzero contributions from the os- 
motic and the elastic terms in the R/R* (or B) expansion 
are smaller by a factor of order (R/R*) 3 . 

Finally, in order to get simply the scaling dependence 
of respectively the osmotic and elastic energy in regimes 2 
and 3 (within unknown numerical coefficients of order 
unity), it is possible to keep only the scaling dependence 
of blobs : £ ~ a<p- 3 / 4 and A ~ N<p-^ 4 /R and to as- 
sign the energy fc^T per blob. Thus, in regime 2, we re- 
cover the very simple estimate of the balance achieved 
by the stretched chains first presented by Alexander jl6| 
and de Gennes p7| . The compressed brush in regime 3 is 
then analogous to a flat brush in presence of a semi-dilute 
solution of monomer volume fraction <Pf, = <P (see Fig- 
ure 2). The parameter R is related to a virtual external 



osmotic pressure 
layer 



9 1 A. 

n ~ &i> which compress the grafted 



i Concluding Remarks 

-In this article, we have given both an SCF and a scal- 
ing study of the structure and the associated free energy 
of polymer brushes in concave geometries. A simple scal- 
ing description can be derived from SCF calculations. The 
essential features are the following : for weak curvatures 
(R > R* ~ 1.3 Lfi a t, where Lfi a t is the thickness of a 
planar brush with the same grafting density), inwardly 
curved brushes are a mere generalization of flat brushes : 
the monomer concentration decreases as we move off from 
the surface. A specific scaling description which would ac- 
count for the geometry and for this slowly decreasing con- 
centration would yield nothing else but constant blobs as 
in the Alexander-de Gennes description of flat brushes. 
This relatively poor scaling description can be related to 
the fact that the chain end distribution spreads over the 
entire layer. However, we argue that a self-similar profile 
of inverted Daoud-Cotton type is not the correct scaling 
description. Even if the variation of the thickness with R 
does not allow to bring to a firm conclusion, this approach 
does not lead to the lowest free energy since not only 
it overestimates the actual free energy but also it yields 
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an increasing monomer concentration which diverges in 
the center for R ~ Lfi at . At large curvatures (keeping 
the grafting density constant), the monomer concentra- 
tion becomes progressively uniform, leading to the struc- 
ture of a squeezed semi-dilute solution where monomer- 
monomer interactions predominate over chain stretching 
energy. This description is consistent with the fact that 
most chain ends are located at the center of the sphere, 
just like in a compressed brush where they are located at 
the edge of the brush. Finally, we also give elements to 
compute numerically the free energy taking into account 
correlations between monomers, by using the SCF results 
for the volume fraction and the chain end distribution. 
We hope that the diagram of the different regimes accord- 
ing to the values of R and the grafting density a will be 
checked experimentally in microemulsions ]49| ]. 
As soon as R ~ Lfi at , the dissymmetry between con- 
vex and concave cases, which originates in the difference 
of available space for the brush, becomes important and 
plays a non negligible role in the stability of micelles and 
polymer decorated membranes. This article brings the es- 
sential elements for the determination of membrane spon- 
taneous curvature, mean and Gaussian rigidity moduli es- 
pecially for decorated vesicles of high curvatures at ther- 
modynamic equilibrium (where the radius is of order of 
100 A [|l2j). This theoretical work should also illuminate 
the issue of the adsorption of polymers on rough surfaces 
consisting of bumps and hollows of high curvatures JlGj ] . 
However, the question of membrane fluctuations with a 
short wavelength A < L fi at still awaits consideration since 
the range of steric interactions between chains is then 
longer than A. The distribution and conformation of grafted 
polymers should then be altered only in the vicinity of the 
grafting points. Indications of this have been given exper- 
imentally by Singh et al. : the influence of surface rough- 
ness on polymer adsorption is smaller when the roughness 
has a very small wavelength lla] . 



The authors wish to thank Alice Nicolas and Bertrand Four- 
cade for useful discussions. We also thank E. Zhulina for point- 
ing out references |24|]2E| ] to us and for raising the argument of 
the stabilizing concentration gradient at the edge of the brush. 

Appendix A: Convolution method for a self- 
consistent field treatment of the concave brush 

We use here an alternate (and faster) method to calculate 
the free energy in the SCF theory. As an example, we treat 
the regime 3. 

We know from the works of Semenov and MWC that the 
potential profile is parabolic, and that it is proportional 
to the volume fraction: 



(r)=#- 



1 r 2 



(39) 



where <P = 3aNa/R is the average volume fraction is cal- 
culated by normalization and where C = 3tt 2 R 2 /8a 2 vN 2 
is found from the "equal-time argument" . In regime 3, the 
volume fraction has some finite value both at the center 
4>{R) = <f> — 9/10(7. Taking advantage of these resu lts, we 
use the convolution method outlined in reference [gOl to 
determine the free energy (in mean- field). From this we 
calculate the two-body excluded volume interaction free 
energy of the brush : 

F s = \ J MR-r) 2 dr^ 2 (r) 



2ttR 3 v f 1 2 
a 3 1 3 



13 
2100 



Defining now (f>fi a t(r) by 

4TrR 2 <j> f i at (r) =4ir(R 



-) a #r) 



(40) 



(41) 



as the density profile obtained by flattening out the con- 
cave brush like a chopped pineapple slice, we compute the 
elastic energy Q : 



FrI, = ±TtR* 



8N 2 a 5 







ir 3 R 5 
2N 2 a 5 



1*. 
10 



dr 



13 

700' 



r 3 d^fi a t(r) 
dr 



C 



(42) 



Gathering all these results, we find the free energy in 
the mean field approach (equation (30)) : 



F = 



2tt R a v 



R 5 



39tt 5 R 7 



20 a 5 N 2 22400 a 7 vN 4 



(43) 



The free-end density profile is then obtained by a simple 
convolution with the volume fraction profile [p0[ : 



g( r o) 



47r(i? - r)' 2 



d(j)fi a 
dr 



i(r) 



r 



Na^/r 



=dr(44) 



where the last factor is the MWC free end distribution for 
a melt brush (step-like volume fraction profile). Through 
integration, we recover equation ((29|). 
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Appendix B: Do lateral fluctuations ruin the 
SCF approach ? 

The SCF approach describes only a projection of the chain 
conformation along the normal to the surface and ignores 
the chain fluctuations parallel to the surface. For a brush 
grafted to a flat interface, the lateral fluctuations can be 
described in very simple terms since the chain performs 
a random walk in these directions : the square-averaged 
position of the n-th monomer from the wall is y = an 1 / 2 
in the case of a melt, and more generally y = em 1 / 2 ^ 1 / 8 
if (j> is the monomer volume fraction. 

Such lateral fluctuations are not hindered in any way 
for a flat brush since the potential is constant in the lateral 
directions. In the case of a concave brush, however, when a 
chain fluctuates laterally, it explores regions located closer 
to the grafting, concave surface, and the potential is some- 
what higher there. In regions where this potential increase 
becomes of order fcgT, we expect a reduction in the lateral 
fluctuations, which in turn represents an additional reduc- 
tion of entropy. In principle, the potential should therefore 
be altered accordingly. We now estimate this effect for the 
regime 3 (where it should be more important) using a 
scaling description of the chain conformations. 

Consider a point at a distance r from the wall where 
the potential has a certain value V(r). Let us now consider 
the potential as a function of the distance z from the first 
point, going perpendicularly to the radius : from geometry, 
the new point is closer to the wall by Ar ~ z 2 /2(R— r). As 
a consequence, the potential is increased by Ar \dV/dr\. 
By considering correlations between monomers, the poten- 
tial equation (H) is modified according to [34}] : V{r') /ksT ■ 
<^/ 4 (r') = (fr^{r'{R)) + 3tt 2 /(8N 2 )(r' 2 (R) - r' 2 ) where 
dr' = ^^(rfdr and 4>{r') = </>(r). Hence, the variation 
of the potential in the perpendicular direction is given by 
(dV/dr')(dr'/dr)Ar, 

In regime 3, the volume fraction is essentially uniform and 
we have r' ~ r^ 1 / 8 thus 

V(r, z) ~ V(r) + fc B T0 1 /4_^!_ (45) 

In the direction of z, the chain evolves in the potential 
given by the above equation. At small length scales, the 
potential increase is weak and does not affect the chain 
conformation. At large length scales, the chain is con- 
fined in this potential. At the crossover between these two 
regimes, the typical distance z that is actually explored 
by the chain is such that the monomers required for a free 
random walk over that distance, each of which being at 
a potential of order V(r, z), correspond to an energy of 
order fc^T altogether. In other words, at a given radial 
position r, the lateral fluctuations of the chain may not 
exceed z given by : 

k B T~ t^l[ V Uz)-V{r)] i.e. z ~ aN 1 ' 2 ^ 1 ^ \?—^- 
a z r 

(46) 

As <p is essentially uniform, the lateral fluctuations depend 
on position only through the last factor [(R — r)/^ 1 / 4 . 



Near the wall, the spontaneous fluctuations are given by : 
z ~ an 1/2 <j)~ 1/a (r < R) (47) 



since this is still lower than z. By contrast, close to the 
center, the lateral fluctuations are given by : 

z~ aA^/^-^l-r/i?] 1 / 4 (R-r<R) (48) 



The energy involved in this effective chain confinement 
is ksT per blob of size z, i.e., per G z ~ (z/a) 2 ^ 1 ' 4 ~ 
N[(R — r)/i?] 1 / 2 monomers. Now, from the SCF chain 
conformation r = Rsm[irn/2N] (if the free end is at the 
center r = R), we deduce that (R-r)/R ~ (N-n) 2 /N 2 . 
Hence, the free energy due to the restriction of the lateral 
fluctuations is given by : 



Fcm ,^r m ^^ nN (49) 

JO 

As a consequence, the free energy that corresponds to the 
reduction of the lateral fluctuations is negligible. The fluc- 
tuations should be included, however, to fully describe the 
chain conformations near the center (for R — r <§; R). 



Appendix C: Cylindrical geometry. 



Following the same approach as in the Introduction, a 
scaling diagram can be established for concave cylindrical 
brushes in terms of the cylinder radius R and the grafting 
density a (Figure 9). It should be noticed that two new 
regimes are present, as compared to the spherical geome- 
try : 



— regime Ibis where the cylinder is so small (i? < aiV 3 / 5 ) 
that mushrooms are confined by the cylinder and stretched 
along the cylinder. This ensures the blob size equal to 
R and their length is ijy ~ N/R 2 / 3 0. 
1/4 regime 4 which is semi-dilute and where chains are not 
laterally compressed. Their length is Ru ~ aiV 1 / 2 ^ 1 / 8 . 



The limit between these two regimes is R ~ l/iV 3 <r 3 . 
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Figure Captions. 

Figure 1: (a) Conformation of a convex brush in the 
Daoud-Cotton scaling description. The geometry ensures 
the equilibrium between elastic and osmotic forces. The 
same picture tr anspos ed to concave surfaces (b) is used 
in references [p"8| , [20|j24|| . The inset shows the monomer vol- 
ume fraction profile. 

Figure 2: Various regimes for a spherical concave grafted 
polymer brush. The inverse of curvature is R and the graft- 
ing density a. A schematic representation of the brush 
is shown using blobs defined in Section 3. This figure is 
only qualitative. (1: mushrooms, 2: weak concave brush, 3: 
compressed brush, 4: collapsed brush, in gray: unphysical 
regions) 



Figure 3: Variations of the free energy per chain F(x) 
for three different values of R/R* for N = 1000, cr = 0.03 
(R*/a = 300.07): from bottom to top R = 2R* (region 2); 
R = R*; R = 0.8R* (region 3). The dashed part of the 
curve corresponds to an unphysical situation [|33[. 



Figure 4: Concentration profile of the brush, (a) Large 
sphere (R = 2 R*, regime 2) : the profile is parabolic and 
the concentration vanishes at some finite distance x* from 
the surface of the sphere (see Eq. (19)). (b) At R = i?*, 
the concentration vanishes precisely at the center of the 
sphere, (c) Small sphere (R = 0.8 R*, regime 3) : trun- 
cated parabolic profile with a finite concentration at the 
center of the sphere (see Eq. (Ed)). 



Figure 5: Variation of the normalized chain end den- 
sity distribution p vs. r$/R for two different radii, R = 
1.27 R* (regime 2) and R = 0.75 R* (regime 3) (N = 1000, 
/ = 600). 



Figure 6: Variation of the solution x* of equation (15) 
vs the normalized curvature R*/R. The inset shows the 
variation of the normalized brush height L/Lfi at vs the 
normalized curvature of the sphere R* /R. 



Figure 7: Variation of the scaled total free energy per 
chain ip (see Eq. (19)) in regimes 2 and 3 vs the normalized 
curvature R*/R. Are also represented the elastic {ipei) and 
osmotic (ipos) parts (where ip = ip os + ip e i) 



Figure 8: Mean square distance between two monomers 
from the same chain (radial component), as a function 
of the number of monomers between them (logarithmic 
scales), (a) In regime 3, for a given average degree of 
stretching \dr/dn\ and volume fraction <fi (thin line), the 
thick line is the scaling picture. At long distances, the de- 
pendence is linear, which reflects the fact that the chain 
is stretched (exponent 1). At smaller distances, the chain 
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takes the conformation it would have in a semi-dilute so- 
lution of volume fraction : it is swollen at small length 
scales (exponent 3/5) and it is Gaussian at length scales 
larger than £ = a</>~ 3 / 4 (exponent 1/2). (b) In regime 4, 
the chains are compressed : the SCF approach does not ap- 
ply and the scaling approach yields the following picture. 
For a given volume fraction cj>, all chains behave alike. 
They are swollen at small length scales (exponent 3/5) 
and Gaussian at length scales larger than £ = a</>~ 3 / 4 (ex- 
ponent 1/2 like in a semi-dilute solution). A G- monomer 
strand spans over the entire interior of the sphere of size 
R and thus larger strands do not extend any further : they 
are confined (exponent 0). 

Figure 9: Various regimes for a cylindrical concave 
grafted polymer brush. The inverse of curvature is R and 
the grafting density a. This figure is only qualitative. (1: 
mushrooms, Ibis: compressed mushrooms, 2: weak con- 
cave brush, 3: compressed brush, 4: semi-dilute brush) 
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